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Abstract 

We present an analysis of particle production yields measured in central Au-Au 
collisions at RHIC in the framework of the statistical thermal model. We demon- 
strate that the model extrapolated from previous analyses at SPS and AGS energy 
is in good agreement with the available experimental data at y^s = 130 GeV im- 
plying a high degree of chemical equilibration. Performing a y 2 fit to the data, the 
range of thermal parameters at chemical freezeout is determined. At present, the 
best agreement of the model and the data is obtained with the baryon chemical 
potential /xg ~ 46 ± 5 MeV and temperature T ~ 174 ± 7 MeV. More ratios, such 
as multistrange baryon to meson, would be required to further constrain the chem- 
ical freezeout conditions. Extrapolating thermal parameters to higher energy, the 
predictions of the model for particle production in Au-Au reactions at yfs = 200 
GeV are also given. 



1 Introduction 



The ultimate goal of ultrarelativistic nucleus-nucleus collisions is to study the properties of 
strongly interacting matter under extreme conditions of high energy density [[I], |2], |3| . Hadron 
multiplicities and their correlations are observables which can provide information on the na- 
ture, composition, and size of the medium from which they are originating. Of particular 
interest is the extent to which the measured particle yields are showing equilibration. The level 
of equilibrium of secondaries in heavy ion collisions can be tested by analyzing the particle 
abundances ||, || ^| [7], [| [Hj or their momentum spectra 0, jl2| . In the first case one estab- 



lishes the chemical composition of the system, while in the second case additional information 
on dynamical evolution and collective flow can be extracted. 

A detailed analysis of the experimental data at SPS energy has shown that hadronic yields 
and their ratios resemble those of a population in chemical equilibrium. Most particle mul- 
tiplicities measured in nucleus-nucleus collisions at the SPS are well consistent with thermal 
model predictions |], [IT], [TTJ . Particle transverse mass distributions, on the other hand, ex- 
hibit slopes showing a structures typical of a thermal source undergoing transverse expansion 



12, [L3| . The slope parameters and the spectra have been well described by the convolution 



of thermal and flow components within the framework of hydrodynamical models Thus 



on the basis of these analyses one concludes that, at the SPS, chemical and thermal equilib- 
rium was indeed achieved at some stage of the collision. The chemical freezeout temperature 
Tf ~ 168 ±2.4 MeV found from a thermal analysis || of experimental data in Pb-Pb collisions 
at SPS is remarkably consistent within error with the critical temperature T c ~ 170 ± 8 MeV 
obtained from Lattice Monte-Carlo simulations of QCD at vanishing baryon density |15 , [16 



Thus, the observed hadrons seem to be originating from a deconfined medium and the chemical 
composition of the system is most likely established during hadronization || [|. The observed 
coincidence of chemical and critical conditions in the QCD medium, if indeed valid, should also 
be seen in heavy ion collisions at higher collision energies, in particular at RHIC. 

In this paper, we present a thermal analysis of particle production measured in Au-Au 
reactions at ^=130 GeV which were reported in [0, ^ |TJ, 0, [E| g|, ^ |2§. We 



will focus on momentum integrated particle yields and their ratios, thereby studying the soft 
component of particle spectra which gives the major contribution to the overall multiplicity. 
We derive the chemical freezeout conditions at RHIC and discuss their uncertainties. We also 
propose an experimental observable which could further constrain the thermal parameters. 
Finally, we provide model predictions for particle production in Au-Au reactions at y/s = 
130 and 200 GeV. 



2 Results 

In the present analysis of RHIC dataQ, we adopt the statistical model which has been used 
previously in the analysis of particle production at SPS and AGS energies ||. This model was 
formulated in the grand canonical ensemble with the baryon number, strangeness and electric 
charge conservation. The partition function contains the contributions from all mesons and 
baryons with masses up to 1.5 and 2.0 GeV, respectively. To account for repulsive interactions 
between hadrons at small distances, an eigenvolume was assigned to all particles. The eigen- 
volume parameter is assumed to be the same for all hadrons and is calculated with a hard core 

1 A summary of data of different collaborations is presented in J27J 
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radius of 0.3 fm. Calculations of relevant observables such as particle densities are performed 
thermodynamically in a self-consistent manner. Thus all these quantities are derived directly 
from the partition function by taking appropriate derivatives. 

Within this framework particle multiplicity ratios depend only on two independent param- 
eters: the temperature T and the baryon chemical potential fiB- All other fugacity parameters 
are determined by the initial conditions requiring strangeness neutrality and electric charge 
conservation. The contribution of resonance decays to the final particle yields of lighter mesons 
and baryons is also included in the analysis. 

The above model is now applied to Au-Au collisions at \fs =130 GeV at RHIC. We use the 
results of the STAR |I7|,|If, 0, 0, 0, PHENIX @§, PHOBOS @, and BRAHMS [2E 



collaborations for different particle multiplicity ratios listed in Table. 1. All these ratios were 
obtained from measured particle yields integrated over transverse momentum. The pseudo- 
rapidity coverage was, however, limited to the central region, —0.5 < rj < 0.5. 

We adjust the free parameters T and fiB to get the best description of the data. The 
calculations are performed under the condition of strangeness neutrality within the 1 unit of 
rapidity covered by the data. Furthermore, we assume that 50 % of all weakly decaying baryons 
are reconstructed in the experiments. For the criterium of the best fit we use the value of x 2 - 
Fig.l shows the contours of constant x 2 i n t ne T — /ib plane. The best fit corresponding 
to the minimum of x 2 — 5.7 for 7 effective degrees of freedom gives T ~ 174 ± 7 MeV and 
fiB — 46 ±5 MeV as the chemical freezeout parameters at this RHIC energy. The corresponding 
strange chemical potential is fi s = 13.6 MeV. Repeating the calculations assuming that as an 
upper limit 100 % of all weakly decaying baryons are reconstructed leads to T ~ 169 ± 7 MeV 
and fis — 48 ± 4 MeV. Since experiments presently do not correct for feed-down from weakly 
decaying baryons, see, e.g., the discussion in |I9| , the reconstruction efficiency will be large, but 
the sensitivity of the calculations to the actual value is small. In the absence of more detailed 
information we use 50 % as reconstruction efficiency. 

The energy per particle of 1.1 GeV obtained with these parameters is consistent with the 
unified freezeout conditions proposed in ||. We note that these parameter values are close to 
those used for a prediction of particle yields at full RHIC energy |28| of T = 168 MeV and 
fi B = 10 MeV.. 

In Fig. 2 we show the comparison of the thermal model results with experimental data. 
One sees in Fig. 2 and also in Table. 1 that the overall agreement is very good. Most of the data 
are reproduced by the model within the experimental errors. The largest deviations are seen 
in the ratios of K*°/h~ and K*°/h~ but still they are on the level of one standard deviation. 
In Table. 1 the predictions of the model for different particle ratios which are still not known 
experimentally are also presented Q As is visible from Fig.l, the minimum in the x 2 surface 
is narrower in fiB direction while more freedom for variation is still possible in T direction 



(see also results presented in [JT0|, |27|, |9], [30]). To further constrain the value of the freezeout 
temperature additional results for particle ratios are required. Of particular relevance would 
be the ratios of multistrange baryons to mesons which are strongly T dependent while being 
only weakly influenced by fiB- In Fig. 3 we show the calculated temperature dependence of the 
f2~/7r~, 2 + //T + and 3~/K~ ratios as possible observables which would constrain further the 
value of T. 

Considering the best values of thermal parameters derived from the present analysis together 



2 These calculations are made assuming a 50% reconstruction efficiency for particles resulting from weak 
baryon decay. This may need to be adapted to specific experimental conditions when comparing ratios involving 
such particles to data. 
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with previous systematics on freezeout conditions in heavy ion collisions in the energy range 



from SIS up to SPS || [31| , one can make predictions for particle production in Au-Au reactions 
at yfs = 200 GeV. Indeed in it was shown that the energy dependence of the baryon chem- 
ical potential can be parameterized phenomenologically as /i^ ~ 1.3GeV(l + yfs J '4.5GeV) _1 . 
This prediction is consistent within error with our present result for RHIC energy. Using this 
parameterization one finds a decrease of baryon chemical potential from fiB — 46 ± 5 MeV 
at y/s = 130 GeV to \i B ~ 29 ± 8 MeV at y/s = 200 GeV. The unified freezeout condition 
of fixed energy/particle ~ 1.1 GeV provides a temperature increase to 177±7 MeV in Au-Au 
reactions at yfs = 200. The freezeout temperature is obviously bounded from above by the 
deconfinement phase transition temperature. The x 2 fit provides only statistical indication of 
the most probable value of thermal parameters. In Table I we show the results of the statistical 
model for different particle ratios at the top RHIC energy. It is clear that the small baryon and 
corresponding strange chemical potentials imply that particle to anti-particle ratios are close 
to unity. Other ratios are also at values near the baryon-free limit. 

Finally, we have to stress that at energies where baryon stopping is dominant, such as at the 
SPS, the comparison of the statistical model with experiment should be done with Att particle 
yields [§, |33|. This procedure strongly reduces the possible influence of dynamical effects on 
particle yields and guarantees that the conservation laws of quantum numbers are fulfilled. In 
the restricted acceptance near midrapidity one needs to account for additional uncertainties in 
the derivation of thermal conditions from the experimental data. It is not excluded that thermal 
parameters at midrapidity could deviate from their values in full phase space as already seen at 
the SPS energy EDI . On the other hand, if the energy is sufficiently high such that data exhibit 
boost-invariant rapidity plateaus, analysis near mid-rapidity should also be little influenced by 
dynamical effects such as, e.g., hydrodynamical flow [34]]. The current energy seems already to 
be close to this regime. 



3 Summary and conclusions 

In conclusion, we have shown that the statistical model in complete equilibrium gives re- 
sults consistent with the experimental data for particle production in Au-Au collisions at 
\/s=130GeV. At this energy the chemical freezeout appears at T = 174±7 MeV and fiB = 46±5 
MeV. The resulting temperature is only slightly higher than that previously found at the SPS 
where for Pb-Pb collisions T = 168 ± 5MeV. This relatively moderate increase of temperature 
could be expected since in the limit of vanishing baryon density the temperature should not 
exceed the critical value required for deconfinement. The substantial decrease of the baryon 
chemical potential from \xb — 270 MeV at SPS to \lb — 45 MeV at RHIC found in these 
calculations shows that at midrapidity we are dealing with a low net baryon density medium. 

We have also discussed possible uncertainties of the results and proposed observables which 
could provide better constraints on thermal parameters. Finally, extrapolating the values of 
thermal parameters to a higher energy, we have made a prediction for particle production in 
Au-Au reactions at ^/i = 200 GeV. 
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Figure 1: Lines of constant x 2 f° r the comparison of experimental particle ratios from RHIC 
experiments and model calculations for a wide range of thermal parameters (T,/i B ). The dot 
represents the parameter set with minimum x 2 located at T = 174 MeV and fi B = 46 MeV. 
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Figure 2: Comparison between RHIC experimental particle ratios and statistical model 
calculations with T = 174 MeV and /xg = 46 MeV. Experimental data are taken from 
171 
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Figure 3: Temperature dependence of different particle ratios calculated with \lb = 50 MeV. 



Ratio Model calculation Expt. Data Expt. Ref. 
y/s = 130 (200) GeV y/s = 130 GeV 



p/p 0.629 (0.752) 0.65 ± 0.07 STAR JTJ 



BRAHMS |2§ 

p/n~ 0.078 (0.089) 0.08 ± 0.01 STAR Jl 



7T-/7T+ 1.007 (1.004) 1.00 ± 0.02 PHOBOS H 

BRAHMS 



K-/K+ 0.894 (0.932) 0.88 ± 0.05 STAR (TJ 

PHENIX || 
PHOBOS 



629 


(0.752) 


0.65 ± 0.07 






n fi/L -i- n D7 






0.60 ± 0.07 






0.64 ± 0.07 


078 


(0.089) 


0.08 ± 0.01 


007 


(1.004) 


1.00 ± 0.02 






0.95 ± 0.06 


894 


(0.932) 


0.88 ± 0.05 






0.78 ± 0.13 






0.91 ± 0.09 






0.89 ± 0.07 


145 


(0.147) 


0.149 ± 0.02 


037 


(0.036) 


0.06 ± 0.017 


032 


(0.033) 


0.058 ± 0.017 


753 


(0.842) 


0.77 ± 0.07 


894 


(0.942) 


0.82 ± 0.08 


040 


(0.039) 




086 


(1.091) 




123 


(0.127) 





PHENIX 22 



PHOBOS |E1 



BRAHMS H 



K-/ir~ 0.145 (0.147) 0.149 ± 0.02 STAR p| 

K*°/h- 0.037 (0.036) 0.06 ± 0.017 STAR f20| 

K^/hT 0.032 (0.033) 0.058 ± 0.017 STAR |20 

A/A 0.753 (0.842) 0.77 ± 0.07 STAR |2| 

S/S 0.894 (0.942) 0.82 ± 0.08 STAR |2T 

A/h- 
A/K*° 
E-/A 

E-/K- 49.8-10- 3 (50.2-10- 3 

S+/A 0.145 (0.142) 

S+/7T+ 6.51-10~ 3 (7.01-10- 3 

0.196 (0.197) 

Q+/Q- 0.898 (0.941) 

n~/ir- 1.47-10" 3 (1.50-10- 3 



Table 1: Comparison of experimental particle ratios and thermal model calculations for T = 174 
MeV, (1b = 46 MeV at ^/s = 130 GeV. Also shown in parentheses are model predictions for 
particle ratios at = 200 GeV with T = 177 MeV, \i B = 29 MeV. 
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